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We  report  the  mechanical  behavior  of  vertically  aligned  carbon  nanotube  films,  grown  on  Si 
substrates  using  atmospheric  pressure  chemical  vapor  deposition,  subjected  to  in  situ  large 
displacement  (up  to  70  pm)  flat-punch  indentations.  We  observed  three  distinct  regimes  in  their 
indentation  stress-strain  curves:  (i)  a  short  elastic  regime,  followed  by  (ii)  a  sudden  instability, 
which  resulted  in  a  substantial  rapid  displacement  burst  manifested  by  an  instantaneous  vertical 
shearing  of  the  material  directly  underneath  the  indenter  tip  by  as  much  as  30  pm,  and 
(iii)  a  positively  sloped  plateau  for  displacements  between  10  and  70  pm.  In  situ  nanomechanical 
indentation  experiments  revealed  that  the  shear  strain  was  accommodated  by  an  array  of  coiled 
carbon  nanotube  “microrollers,”  providing  a  low-friction  path  for  the  vertical  displacement. 
Mechanical  response  and  concurrent  deformation  morphologies  are  discussed  in  the  foam-like 
deformation  framework  with  a  particular  emphasis  on  boundary  conditions. 


I.  INTRODUCTION 

Considerable  efforts  have  been  dedicated  to  explore 
the  deformation  mechanisms  of  vertically  aligned  carbon 
nanotube  (VACNT)  forests,  in  part  motivated  by  their 
wide  range  of  potential  applications  in  areas  such  as 
energy  dissipation  devices,  electrical  interconnects,  ther¬ 
mal  interface  materials,  microelectromechanical  systems, 
and  microelectronics.'^®  To  date,  instrumented  indenta¬ 
tion  has  been  a  common  method  for  testing  the  mechanical 
behavior  of  VACNT  films  on  substrates.  However,  due  to 
the  limitations  of  the  displacement  actuators,  most  of  the 
existing  studies  are  limited  to  shallow  indentation  depths 
(and  strains),  varying  from  a  few  hundred  nanometers  to 
several  micrometers.  Additionally,  while  VACNT  forest 
films  may  have  millimeter-sized  lateral  and  vertical 
dimensions,  they  are  composed  of  individual  nanotubes 
with  diameters  in  the  nanometer  range,  which  drives 
their  mechanical  response  to  be  distinct  from  monolithic 
materials.  Hence,  the  accurate  estimation  of  the  contact 
area  between  the  VACNT  film  and  the  commonly  used 
parabolic^’*  and  pyramidal^  indenter  tip  geometries. 
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necessary  for  indentation  data  analysis,  poses  a  significant 
challenge. 

In  this  work,  we  utilize  flat-punch  diamond  indenter 
tips,  resulting  in  a  constant  contact  area  with  the  sample, 
to  conduct  large  displacement  indentations  into  VACNT 
films  grown  using  atmospheric  pressure  chemical  vapor 
deposition  (APCVD).  Using  both  on-edge  in  situ  and  in- 
bulk  ex  situ  flat-punch  indentations,  we  examine  the 
combined  quantitative  information  (load-displacement) 
along  with  the  concomitant  morphological  changes.  To 
explore  the  effect  of  lateral  boundary  constraints,  these 
findings  are  compared  with  our  previous  report^  on  uni¬ 
axial  compressions  of  the  same  APCVD-grown  VACNTs 
patterned  into  cylindrical  bundles. 

Mechanical  response  of  VACNTs  is  a  complex  phe¬ 
nomenon  occurring  across  multiple  length  scales  due  to 
its  hierarchical  structure.  The  collective  behavior  of  these 
materials  is  expected  to  rely  heavily  on  the  properties  of 
the  individual  carbon  nanotubes  (CNTs),  as  well  as  on 
the  variations  in  the  collective  intertube  interactions  and 
inherent  property  gradients  of  the  microstructure,  which 
in  turn  are  dictated  by  their  synthesis  techniques. 

A  wide  range  of  mechanical  behavior  has  been  reported  in 
VACNTs  as  a  result  of  such  variations  in  their  micro¬ 
structure.  These  include  modulus  and  buckling  strengths 
ranging  from  sub-MPa'^’'"'  to  tens  of  MPa'®^'^  to  GPa,*’'* 
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as  well  as  a  complete  gamut  of  recoverability — from  nearly 
100%  recovery  even  after  large  deformations^’^’ to 

011  1701  00 

permanent  distortion  even  at  modest  strains.  >  ’  ’  ’ 
Such  drastic  differences  in  properties  make  it  impractical 
to  compare  VACNTs  synthesized  by  nonidentical  growth 
techniques.  In  this  work  we  focus  on  VACNTs  grown  by 
APCVD  technique^^  and  describe  their  response  under 
indentation.  These  results  are  contrasted  with  our  pre¬ 
vious  work  on  uniaxial  compressions  of  30  pm  diameter 
x30  pm  tall  cylindrical  bundles.®  The  APCVD  VACNT 
bundles  are  known  to  exhibit  three  distinct  regimes  in  their 
compressive  stress-strain  curves:  (i)  a  short  initial  elastic 
section,  followed  by  (ii)  a  sloped  plateau  with  character¬ 
istic  wavy  features  corresponding  to  localized  buckle 
formation,  and  (iii)  a  densification  characterized  by  a  rapid 
stress  increase — all  of  which  liken  their  deformation  to 
that  of  a  typical  foam.^’®’"  Unlike  foams,  however,  the 
compressive  strain  in  VACNTs  is  accommodated  entirely 
via  the  formation  of  localized  folds  or  buckles  forming 
progressively  along  their  lengths.  •  ■  •  •  •  During 
compression  of  our  particular  APCVD  VACNT  system, 
such  localized  buckles  fully  recovered  upon  load  release 
even  from  relatively  large  strains  of  60-80%  (i.e.,  the 
onset  of  densification)  under  fast  displacement  rates  of 
1000  nm/s.  In  addition,  repeated  cycling  also  revealed 
a  distinctive  hysteresis  behavior  for  the  VACNT  bundle  in 
compression,  suggesting  energy  dissipation  in  every 
cycle.®  However,  since  these  properties  were  mainly 
documented  for  short  height  APCVD  VACNTs  under 
uniaxial  compression,  it  is  still  largely  unknown  whether 
they  demonstrate  similar  lucrative  mechanical  properties 
under  different  loading  and  boundary  conditions. 

In  this  work,  we  focus  on  the  deformation  response  of 
I42-pm-tall  APCVD-grown  VACNT  films  under  large 
displacement  flat-punch  indentations  and  compare  it  with 
that  of  uniaxial  compressions  of  —30  x  30  pm  (height  x 
diameter)  cylindrical  pillars  patterned  into  the  same 
starting  material.  Thus,  we  directly  compare  the  influence 
of  the  boundary  conditions  posed  by  the  presence  or 
absence  of  the  external  matrix.  We  use  both  an  in  situ 
nanoindentation  methodology  to  observe  the  on-edge 
deformation  in  real  time  using  a  custom-built  in  situ 
nanomechanical  deformation  instrument,  SEMentor,^®  as 
well  as  more  traditional  ex  situ  indentation  methods. 
Recent  work  by  our  group®’' and  others showed 
the  efficacy  of  conducting  in  situ  compressions  inside  the 
scanning  electron  microscope  (SEM),  which  allows  for 
uninterrupted  observation  of  the  real-time  evolution  of 
deformation  morphology  while  simultaneously  recording 
load  versus  displacement  data.  Unlike  compression 
experiments,  where  the  indenter  tip  applies  an  axial  load 
to  an  individual  structure  readily  visible  in  SEM,  in  situ 
indentation  experiments  have  to  be  conducted  on  the 
edge  of  the  sample,  as  opposed  to  in  bulk,  to  facilitate 
uninhibited  inspection  of  the  material  cross  section. 


Although  boundary  conditions  and  constraints  associated 
with  such  on-edge  tests  are  different  from  those  during 
in-bulk  indentations,  they  are  instrumental  in  providing 
useful  information  on  morphological  evolution  in  deform¬ 
ing  VACNTs,  not  easily  obtainable  by  other  methods. 

II.  MATERIALS  AND  METHODS 

A.  CNT  growth 

VACNT  films  with  thicknesses  between  25  and  142  pm 
were  grown  using  APCVD  technique  in  a  commercial 
chemical  vapor  deposition  (CVD)  system  (Black  Magic 
Pro  4";  Aixtron  SE,  Herzogenrath,  Germany)  as  described 
in  detail  in  Ref.  23.  A  trilayer  catalyst  of  Ti  (30  nm)/Al 
(10  nm)/Ee  (3  nm)  was  deposited  by  electron  beam 
evaporation  on  Si  substrates  to  seed  the  growth  of  the 
VACNT  films.  A  gas  mixture  of  160/100/7500  seem  of 
C2H2/H2/N2  was  used  to  maintain  the  chamber  pressure  at 
750  mbar  during  growth,  which  occurred  at  —750  °C. 
Multiwall  CNTs,  with  an  average  diameter  of  8.8  ±  2. 1  nm 
(average  ±  standard  deviation)  and  average  mass  density 
(CNT  mass  divided  by  volume)  of  —80  mg/cm^,  were 
produced  by  this  process. 

B.  Ex  situ  indentation  tests 

Two  different  VACNT  films  with  thicknesses  of  —142 
and  25  pm  were  selected  for  large  displacement  cyclic 
indentation  tests.  The  indentation  experiments  were  per¬ 
formed  using  the  XP  module  of  Agilent’s  nanoindenter 
G200  (Santa  Clara,  CA)  with  adjustable  software  control 
methods  as  described  in  Ref.  1 1 .  Tests  were  performed  in 
air  using  a  custom-made  cylindrical  diamond  fiat  punch 
with  —120  pm  diameter  and  —80  pm  height.  Eor  the 
142-pm-thick  VACNT  sample,  indentations  were  performed 
under  a  constant  displacement  control  to  varying  indentation 
depths  up  to  a  maximum  penetration  depth  of  around  70  pm 
(restricted  by  the  height  limitations  of  the  diamond  fiat 
punch).  Eor  the  thinner  samples,  the  maximum  depth  was 
limited  to  90%  of  the  film  thickness.  Tests  were  conducted  in 
the  interior  of  the  samples  (i.e.,  away  from  the  sample  edge) 
at  three  different  constant  displacement  rates:  10,  100,  and 
1000  nm/s.  Typically,  five  unload-reload  cycles  were 
performed  at  each  displacement  level.  No  hold  time  was 
applied  at  the  maximum  loads  to  maintain  the  cyclic  nature  of 
the  tests.  Minimum  of  10  tests  were  conducted  at  each 
displacement  rate.  The  indents  were  spaced  at  least  500  pm 
apart  to  eliminate  any  possible  proximity  effects. 

C.  In  situ  SEM  indentation  experiments 

To  visualize  the  morphological  changes  in  the  VACNT 
samples  during  indentation  and  correlate  them  with  the 
concurrent  load-displacement  data,  in  situ  tests  were 
conducted  in  a  custom-built  mechanical  deformation  in¬ 
strument,  SEMentor,^^  composed  of  a  nanomechanical 
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dynamic  contact  module  (Agilent  Corp.,  Santa  Clara,  CA) 
inside  a  SEM  (Quanta  200;  FEI,  Hillsboro,  OR).  Tests  in 
SEMentor  were  conducted  with  a  conductive  diamond  flat 
punch  with  a  rectangular  flat  cross  section  of  —60  x  80 
pm.  The  loading  axis  in  SEMentor  is  inclined  at  —86° 
with  respect  to  the  electron  beam,  thus  allowing  contin¬ 
uous  observation  of  the  deformation  morphology  of  the 
VACNT  film  cross  section  during  the  on-edge  in  situ 
experiments.  Importantly,  the  outer  constraints  in  on-edge 
indentations  are  different  from  those  during  in-bulk  tests, 
and  therefore,  their  stress-strain  behavior  is  expected  to  be 
different,  as  well.  SEM  observations  were  recorded  as 
a  video  file  at  30  frames/s  capture  rate  and  synchronized 
with  the  indentation  data  to  provide  a  real-time  correlation 
between  each  video  frame  and  the  corresponding  position 
on  the  load-displacement  curve  (see  video  files  S 1  and  S2, 
shown  at  15  and  10  times  their  original  speeds,  respec¬ 
tively,  in  the  supporting  online  material).  SEMentor  tests 
were  conducted  to  a  maximum  penetration  depth  of  30  pm 
(instrument  limit).  In  addition,  while  nominally  the  experi¬ 
ments  performed  in  both  G200  and  SEMentor  are  identi¬ 
cal,  one  notable  difference  is  that  SEMentor  tests  are 
conducted  in  a  vacuum  environment  while  those  in  G200 
are  conducted  in  air.  Eurther,  in  SEMentor  experiments, 
the  samples  are  constantly  exposed  to  the  electron  beam, 
and  they  are  oriented  horizontally  such  that  gravity  is  acting 
perpendicularly  to  the  compression  axis.  Due  to  these  subtle 
differences  between  the  two  instmments,  the  SEMentor- 
generated  results  were  used  only  for  visualization  purposes 
in  this  study,  while  all  data  analyses  were  performed  on  tests 
conducted  in  air  in  the  G200  nanoindenter. 

D.  Data  analysis 

The  applied  load,  P,  and  measured  displacement,  h, 
were  corrected  for  machine  compliance  following  the 
procedure  outlined  in  detail  in  Ref.  11.  We  assumed 
negligible  friction  between  the  indenter  sidewalls  and 
the  VACNT  matrix — a  reasonable  approximation  based 
on  previous  reports.  ’  The  initial  loading  portion  in 
a  typical  flat-punch  indentation  experiment  often  suffers 
from  misalignment  issues  between  the  indenter  tip  and 
the  sample  surface  requiring  special  care  in  reliable  extrac¬ 
tion  of  useful  properties.  On  the  other  hand,  applying  the 
Hertzian  contact  model  to  the  initial  unloading  segment 
of  the  measured  load-displacement  curve  permits  a  rea¬ 
sonable  estimation  of  the  Young’s  modulus^®’^°: 

1  _  1  -  I -vj 

Eeff  Es  Ei 

(1) 

where  Eeff  denotes  the  effective  modulus  of  the  combined 
indenter-specimen  system,  S(=dP/dh)  is  the  stiffness 
measured  from  the  slope  of  the  initial  30%  of  the 


S  _  S 
2  VA  2a 


unloading  load-displacement  curve,  v  and  E  are  the 
Poisson’s  ratio  and  the  Young’s  modulus,  respectively, 
and  the  subscripts  s  and  i  refer  to  the  specimen  and  the 
indenter,  respectively.  Eollowing  the  common  nanoinden¬ 
tation  protocols,  we  chose  the  modulus  of  1041  GPa  and 
Poisson’s  ratio  of  0.07  for  the  diamond  indenter  tip.  Since 
the  APCVD-grown  VACNT  bundles  have  been  shown  to 
exhibit  foam-like  deformation  with  minimal  lateral  expan¬ 
sion  under  compressive  load,^  we  assume  a  vanishing 
Poisson’s  ratio  of  v  =  0  for  the  VACNTs  as  has  also  been 
suggested  previously.^ ^  We  note  that  this  choice  may  be 
an  approximation  since  the  deformation  in  VACNTs  is 
highly  localized.  However,  the  error  introduced  is  not 
expected  to  be  significant  since  the  accuracy  of  depends 
only  weakly  on  v  [see  Eq.  (1)];  e.g.,  a  change  in  v  from 
0  to  0.3  causes  only  a  5.3%  error  in  E^. 

We  also  note  that  Eq.  (1)  treats  the  overall  material 
behavior  as  that  of  an  isotropic  continuum  medium, 
according  to  the  inherent  assumptions  of  the  Hertz’s  theory. 
VACNTs,  on  the  other  hand,  demonstrate  varying  degrees 
of  anisotropy  at  each  level  of  its  hierarchical  microstmcture 
(see  Eig.  1),  rendering  our  calculations  of  indentation 
moduli  an  approximation  rather  than  an  exact  value. 
Nevertheless,  isotropic  continuum  framework  has  been 
previously  utilized  in  developing  the  constitutive  relations 
in  VACNTs"  and  foams and  appears  to  have  accurately 
captured  the  qualitative  features  of  outer  deformation  pro¬ 
file  and  the  stress-strain  curves.  Despite  the  discrete 
constituents  and  apparent  vertical  alignment  at  low  magni¬ 
fication,  the  continuum  foundation  was  motivated  by  the 
nearly  isotropic  network  of  CNTs  as  revealed  by  images  at 
or  above  magnifications  of  240,000x  [Eig.  1(c)]  and  above. 

We  define  the  percentage  recovery  (R)  as  the  displace¬ 
ment  recovered  at  the  end  of  each  cycle  with  respect  to  the 
maximum  displacement,  i.e.,  R  =  where  /j^ax  is 

“max 

the  maximum  displacement  at  the  end  of  loading  and 
/lunioad  IS  the  displacement  after  unloading  to  10%  of  the 
maximum  load  in  each  cycle. 

The  loss  coefficient,  q  (a  dimensionless  quantity), 
measures  the  degree  to  which  a  material  dissipates  energy 
and  is  calculated  as^^: 


= 


AUi 

iTlUi 


U  —  /  ^indd^ind 


1  af 


2  E. 


AU  (p  ; 


(2) 


where  U i  is  the  elastic  energy  stored  in  the  material  when 
it  is  loaded  elastically  to  a  stress  dmax  in  the  1st  cycle  and 
AUi  is  the  energy  dissipated  in  the  ith  load-unload  cycle 
[see  Eig.  5(d)]. 

To  compare  the  indentation  response  of  VACNTs  to  its 
stress-strain  behavior  under  compression,  we  define  the 
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Indent  on 
sample  bulk 


FIG.  1.  SEM  images  reveal  the  hierarchical  moiphology  of  the 
(a)  141.5-pm-thick  VACNT  films  (magnification  260x),  which  consist 
of  (b)  nominally  vertical  aligned  CNTs  visible  at  a  lower  magnification 
of  30kx,  and  (c)  a  complex  intertwined  network  seen  at  higher 
magnifications  of  240kx.  SEM  pictures  are  taken  at  a  60°  tilt  angle, 
(c  inset)  Individual  multi  walled  CNTs  of  outer  diameter  8.8  ±  2.1  nm 
(average  ±  standard  deviation)  are  visible  in  the  TEM  image. 


average  indentation  stress  (ainu)  and  strain  (Sind)  under  a 
flat-punch  cylindrical  indenter  of  contact  diameter  2a  and 
contact  area  A  as  given  below. In  Eq.  (3),  the  indentation 
response  is  assumed  to  free  of  any  substrate  effects: 


P 

a.„d-- 


P 

71(3^ 


_  h 

^ind 

2a 


(3) 


E.  Viscoelastic  characterization 

The  viscoelastic  properties  of  the  VACNT  film  were 
measured  following  the  procedure  outlined  in  Ref.  1 1 .  In 
this  method,  the  indenter  is  loaded  into  the  sample  in  air  at 
a  constant  displacement  rate  of  100  nm/s  up  to  a  specified 
indentation  depth  where  the  indenter  head  is  oscillated  at 
~8  nm  amplitude  across  a  range  of  frequencies  from  1  to 
50  Hz.  The  cutoff  frequency  of  50  Hz  is  dictated  by  the 
instmmental  limit,  as  detailed  in  Refs.  11  and  34.  The 
procedure  was  repeated  at  seven  different  constant  in¬ 
dentation  depths:  from  5  to  70  pm. 

Viscoelastic  materials  are  commonly  characterized  by 
their  storage  (£')  and  loss  (£")  moduli,  where  the  former 
represents  the  stored  energy  or  the  elastic  response  and  the 
latter  corresponds  to  the  amount  of  energy  dissipated,  or 
the  viscous  response,  as  well  as  their  ratio — tan5.  Assum¬ 
ing  linear  viscoelastic  behavior,  these  terms  can  be 
computed  following  the  calculations  described  in  Refs. 
34-37  as  follows: 


E'  =  k' 

2|3  ’ 

E"  =  k" 


2P  VA 


tan5  =  — 


Mo 

Mo 


COS  (p  — 

sincp  — 


Mo 

Mo 


COS  (pair. 


smcpa 


(4) 


Here,  k'  and  k"  are  the  storage  and  loss  stiffnesses  of  the 
sample,  obtained  by  finding  the  real  and  complex  parts, 
respectively,  of  the  stiffness  differences  between  oscillat¬ 
ing  the  indenter  head  on  the  sample  at  a  fixed  displacement 
and  in  air  at  the  same  raw  displacement,  and  P  is  a  constant 
(  =  1  for  a  flat-punch  indenter).  Fq  and  mq  are  the  load  and 
displacement  oscillation  amplitudes,  respectively,  and  cp  is 
the  phase  angle  between  the  load  and  displacement 
oscillations.  We  note  that  the  accuracy  in  the  values  of 
E'  and  E'  in  Eq.  (4)  can  be  affected  by  several  factors: 
uncertainties  in  Poisson’s  ratio,  since  it  may  be  frequency 
dependent,  and  ambiguity  in  contact  area,  especially  at 
shallower  indentation  depths,  where  full  contact  may  not 
have  been  established.  On  the  other  hand,  calculation  of 
tan5  is  independent  of  the  contact  area  and  thus  is  ideally 
suited  as  a  measure  of  the  viscoelasticity  of  the  indented 
material.*’^^ 
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III.  RESULTS  AND  DISCUSSION 
A.  VACNT  microstructure 

Figure  1  shows  the  SEM  and  transmission  electron 
microscopy  (TEM)  images  at  different  magnifications  of 
a  ~142-pm-thick  VACNT  film  on  a  Si  substrate  grown 
using  the  APCVD  synthesis  technique  (see  Sec.  II).  All 
SEM  images  [Eigs.  l(a)-l(c)]  were  taken  at  a  60°  tilt  to 
visualize  both  top  and  edge  surfaces  of  the  sample.  These 
figures  show  the  complex  hierarchical  nature  of  the 
VACNT  microstructure  with  their  distinct  organizational 
details  across  multiple  length  scales.  Thus,  while  the 
VACNT  appears  as  a  continuous  film  at  lower  magnifica¬ 
tions  of  260x  [Eig.  1(a)],  the  nominally  vertical  alignment 
of  CNTs  growing  perpendicularly  to  the  support  substrate 
in  the  VACNT  bundle  becomes  visible  at  a  higher 
magnification  of  30,000x  [Eig.  1(b)].  Increased  magnifi¬ 
cation  of  240,000x  reveals  significant  intertwining  in  the 
long  curved  lengths  of  individual  CNTs,  revealing  the 
isotropic  CNT  network. As  a  result,  in  their  as- 
grown  states,  certain  segments  of  the  CNT  forest  appear  to 
be  prebuckled/prebent,  and  the  favorable  contact  energy 
between  the  tubes  (van  der  Waals)  is  thought  to  balance 
the  bending  strain  energy  of  their  arrangement,  resulting  in 
a  stable  low  energy  configuration.^^  The  diameters  of  the 
individual  multiwalled  (3-6  walls)  CNTs  (8.8  ±2.1  nm) 
are  visible  at  the  nanometer  length  scale  shown  in  the  TEM 
image  in  Eig.  1(c)  inset.  Such  hierarchy  in  the  VACNT 
microstructure  governs  its  overall  mechanical  behavior, 
which  is  a  result  of  not  only  the  properties  of  individual 
CNTs  but  also  of  their  complex  mutual  interactions  and 
distribution  throughout  the  array. 

Eigure  1(a)  shows  the  residual  imprints  of  indentations 
(indicated  by  arrows):  one  ex  situ,  performed  in  the 
interior  of  the  sample  away  from  the  edges,  and  two  in 
situ  indentations,  performed  on  a  sample  edge  inside  the 
SEMentor  using  a  rectangular  —60  x  80  pm  fiat  punch 
(see  Sec.  II).  As  can  be  seen  clearly  from  this  image,  the 
indentations  left  a  significant  imprint  on  the  surface  of  the 
VACNT  film,  suggesting  marginal  postmortem  elastic 
recovery.  This  is  in  contrast  with  the  nearly  full  recover¬ 
ability  exhibited  by  the  same  VACNTs  patterned  into 
micrometer-sized  bundles  and  compressed  to  the  same 
nominal  strain.^  This  permanent  shape  change  underneath 
the  indenter  appears  to  have  occurred  along  the  vertical 
edges  via  a  single  shear  offset  (see  also  Eig.  3).  Notably, 
although  the  two  indents  on  the  sample  edge  shown  in 
Eig.  1(a)  are  only  60  pm  apart  from  one  another,  we  see  no 
evidence  of  any  infiuence  of  the  earlier  (left)  indent  on  the 
latter  (right).  This  observation  suggests  that  the  lateral 
extent  of  the  plastic  zone  size  in  VACNTs  under  the 
indenter  is  much  narrower  as  compared  with  the  rounded 
continuum  indentation  stress  field  contours  spanning  the 
regions  as  far  as  one  indenter  diameter  away  from  the  fiat- 
punch  indenter  tip  in  monolithic  materials. This  further 


underlines  the  importance  of  understanding  the  hierar¬ 
chical  nature  of  VACNT  materials  as  their  postindent 
morphological  signature  is  clearly  of  nonmonolithic  nature, 
while  its  stress-strain  response  is  akin  to  that  of  an  open¬ 
cell  foam. 

B.  In  situ  on-edge  flat-punch  indentations 

To  correlate  the  morphological  changes  underneath 
the  indenter  with  the  details  of  the  indentation  load- 
displacement  curve,  a  series  of  in  situ  indentations  were 
performed  in  SEMentor  using  the  procedures  described 
in  Sec.  II.  Eigure  2(a)  shows  two  representative  load- 
displacement  curves,  with  the  respective  indentation  stress 
values  plotted  on  the  secondary  y  axis.  These  correspond 
to  two  of  the  three  indentations  shown  in  Eig.  1 :  one  on  the 
edge  and  another  in  the  sample  interior.  Both  tests  were 
performed  at  a  nominal  displacement  rate  of  100  nm/s. 
This  plot  indicates  that  while  the  curves  are  similar  in 
shape,  the  on-edge  regions  deform  at  significantly  lower 
loads  than  in-bulk  locations  at  equivalent  indentation 
depths.  This  result  is  not  surprising  due  to  the  less  restrictive 
boundary  conditions  in  the  on-edge  setup.  In  both  cases,  the 
load  first  increases  steeply  up  to  an  instability  manifested  by 
a  large  displacement  burst  of  —20  pm,  after  which  the 
indenter  tip  is  unloaded  [Eig.  2(a)].  Video  frames  from  the 
indent-on-edge  setup  are  provided  in  supporting  online 
video  (file  SI),  with  several  representative  snapshots  shown 
in  Eigs.  2(b)-2(d).  As  revealed  by  these  images,  the  first 
buckle  marked  by  double  arrows  in  Eig.  2(c)  forms  close  to 
the  substrate  during  the  initial  loading  at  the  load  of  —6  mN. 
This  buckle  is  initially  observed  only  in  the  region  just 
underneath  the  indenter  tip — the  lateral  features  outside  of 
this  region  marked  with  single  arrows  appear  after  the 
displacement  burst  (see  supporting  online  video  SI).  In¬ 
terestingly,  the  first  buckle  location  was  always  —9  pm 
above  the  substrate.  Other  tests  on  VACNT  samples  of  dif¬ 
ferent  film  thicknesses  have  also  shown  a  similar  response 
where  the  first  buckle  invariably  formed  close  to  the  hot- 
tom  substrate.  This  indicates  that  the  presence  of  a  stiff 
substrate  results  in  localized  buckle  formation  under  both 
compression^’^’' and  nanoindentation. 

Upon  further  loading  the  buckle  extends  laterally 
[region  marked  with  the  single  arrows  in  Eig.  2(c)]  until 
the  load  of  —7.3  mN,  where  a  sudden  and  extensive 
displacement  burst  of  —20  pm  occurs.  SEM  images 
obtained  immediately  after  this  burst  (video  capture  frame 
rate  was  not  fast  enough  to  catch  the  details  during  the 
instantaneous  burst)  reveal  that  a  portion  of  VACNT  film 
sheared  off  nearly  perfectly  vertically  along  the  edges  of 
the  indenter  tip  [Eig.  2(c)],  and  new  buckle  lines  are  visible 
next  to  the  initial  buckles. 

Closer  inspection  of  the  sheared-off  region  conveys 
several  noteworthy  features.  Eirst,  the  deformation  appears 
to  be  highly  localized,  occurring  only  at  the  shear  offset 
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FIG.  2.  (a)  Comparative  indentation  load  (and  stress)-displaeement  curves  at  100  nm/s  loading  rate  of  in  situ  SEM  indentation  tests  conducted  on  the 
sample  edge  versus  in  the  interior  of  the  sample  using  a  60  x  80  pm  rectangular  flat-punch  indenter.  (b,  c,  and  d)  For  the  indent  on  edge,  the  indentation 
process  was  recorded  as  a  video  file  from  which  individual  picture  frames  were  extracted  corresponding  to  a  particular  event  during  the  loading¬ 
unloading  cycle  (as  indicated  on  the  load-displacement  graph).  The  double  arrows  in  (c)  mark  the  initiation  of  the  first  buckle  at  a  load  of  around  6  mN. 
The  wide  lateral  buckles  marked  by  single  arrows  are  formed  after  the  first  buckle  and  increase  in  size  with  loading,  (d)  Upon  unloading,  these  lateral 
buckles  disappear,  (d  inset  top)  The  large  displacement  burst  results  in  a  vertical  shear  of  the  VACNT  structure  along  the  indenter  edges.  Note  the  clean 
edge  of  the  shear  wall,  (d  bottom  insets)  The  shear  appears  to  have  been  carried  by  a  series  of  CNT  “microrollers,”  which  act  as  effective  lubricants 
protecting  the  rest  of  the  VACNT  structure  from  further  damage.  All  SEM  images  except  the  bottom  two  insets  for  (d)  were  taken  at  86°  tilt  angle.  The 
bottom  two  insets  for  (d)  were  at  60°  tilt  angle. 


regions  and  the  buckles  formed  at  the  bottom  of  the  sample. 
The  remaining  portions  of  the  VACNT  film  appear  com¬ 
pletely  unaffected  [Fig.  2(d)  top  inset].  Second,  the  shear 
appears  to  have  been  carried  by  a  series  of  tangled 
CNTs,  or  “microrollers,”  seen  in  the  magnified  inset 
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images  in  Fig.  2(d),  whereby  originally  vertical  strands 
of  several  entangled  CNT  collectively  coiled  themselves 
into  helical  structures.  Three  pairs  of  such  coiled  strands 
of  CNTs  can  be  seen  on  either  side  of  the  sheared  region. 
We  find  the  length  of  the  coiled  CNTs  in  each  microroller 
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to  be  ~6  pim,  and  thus,  the  total  displacement  carried  by 
the  rollers  is  —18  pm.  This  distance  is  in  reasonable 
agreement  with  the  displacement  burst  extent  of  —20  pm 
shown  in  Fig.  2(a).  We  hypothesize  that  these  microrollers 
provide  a  low-friction  path  during  the  shear  process. 

The  shear  event  and  the  formation  of  the  CNT  micro¬ 
rollers  [Figs.  2(c)  and  2(d)]  are  unique  characteristics  seen 
only  in  the  intertwined  VACNT  systems.  Neither 
foams,"^°’"^'  where  the  deformation  is  confined  only  to 
regions  directly  beneath  the  indenter,  nor  noninteracting 
VACNTs,'*^  where  such  shear  would  propagate  without 
any  microroller  formation,  show  the  unique  features  seen 
in  Fig.  2.  As  shown  before  in  Fig.  1,  the  hierarchical  design 
of  these  materials  results  in  considerable  interlinking 
between  the  individual  struts.  Thus,  in  spite  of  their 
nominally  vertical  structure,  the  CNTs  along  the  shear 
plane  of  the  VACNT  film  can  perform  as  highly  efficient 
lubricants  by  collectively  coiling  themselves  and  thereby 
carrying  the  shear  strain  in  response  to  axial  loading.  The 
effectiveness  of  this  “lubrication  process”  is  evidenced  by 
the  very  rapid  nature  of  the  burst  and  the  resulting  shear 
event  under  the  indenter.  Comparing  SEM  images  at 
maximum  load  and  after  unload  [Figs.  2(c)  and  2(d)] 
shows  that  this  shear  event  caused  permanent  damage  to 
the  material  structure,  as  unloading  the  indenter  tip  results 
in  minimal  recovery.  Upon  unloading,  the  CNT  micro¬ 
rollers  did  not  uncoil,  and  the  only  observed  difference  in 
the  morphology  of  the  before  and  after  deformation 
material  was  that  the  buckle  contours,  which  extended 
laterally  beyond  the  sheared  block,  disappeared  upon 
unloading. 

C.  Ex  situ  in-bulk  large  displacement  indentations 

To  analyze  the  deformation  beyond  the  initial  large 
burst,  ex  situ  indentation  tests  were  conducted  in  air  in 
the  Agilent  G200  nanoindenter  (see  Sec.  II)  to  larger 
depths  of  —70  pm  in  the  interior  of  the  as-grown  VACNT 
film  away  from  fhe  edges.  This  maximum  indenfafion 
depfh  was  chosen  based  on  fhe  heighf  of  fhe  cylindrical 
punch  (—80  pm)  and  to  avoid  fhe  influence  of  fhe  rigid  Si 
subsfrafe.  Figure  3(a)  shows  a  fypical  indenfafion  load 
(slress)-displacemenl  curve  for  a  represenfafive  fesf  con¬ 
ducted  al  fhe  displacemenl  rate  of  10  nm/s  up  to  a  final 
displacemenl  of  70  pm  on  fhe  142-pm-lhick  VACNT  film. 
Similar  to  fhe  compression  experimenls  on  fhe  same 
VACNT  system,^  three  distinct  regions  are  apparent:  (i)  a 
short  elastic  regime,  followed  by  a  (ii)  sudden  instability, 
which  results  in  a  40%  load  drop  at  — 4.3-pm  displacement 
and  a  subsequent  (iii)  sloped  plateau  region  with  a  small 
positive  slope  starting  at  the  displacement  of  —11  pm, 
where  the  load  increases  more  gradually  with  displacement 
and  contains  several  undulations.  We  calculate  at  the 
point  of  instability  to  be  —1.6  MPa,  a  value  nearly  40% 
below  the  buckling  stress  of  —2.5  MPa  exhibited  by  this 
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FIG.  3.  (a)  Five  load-unload  cycles  during  indentation  (in  bulk)  on 
a  141.5-|J.m-thick  VACNT  film  at  10  nm/s  loading  rate  using  a  120-)rm- 
diameter  flat-punch  indenter  showing  three  distinct  regimes;  elastic, 
drop  in  load  at  a  critical  instability  stress  ains,  and  plateau  regime, 
(a  inset)  Successive  load-unload  cycles  to  the  same  maximum  displace¬ 
ment  reveal  a  gradual  drop  in  recovery  in  every  cycle,  suggesting 
progressive  damage  accumulation  in  the  material.  The  corresponding 
indentation  stress  values  are  shown  along  the  secondary  y  axis,  (b) 
Indentations  at  different  displacement  rates  showing  a  large  displace¬ 
ment  burst  and  more  pronounced  buckling  signatures  at  the  faster  rates. 
The  inset  images  show  the  residual  impressions  of  two  different 
indentations  conducted  at  a  rate  of  100  nm/s,  one  when  the  indenter  is 
unloaded  just  after  the  burst  versus  when  it  is  unloaded  from  a  maximum 
displacement  of  —65  pm. 


material  system  when  patterned  into  30  x  30  pm,  heighf  x 
diameter,  bundles  under  compression  in  fhe  same  inslm- 
menl.^  Probable  causes  for  Ihis  difference  belween  the 
indentation  and  compression  tests  are  detailed  in  fhe 
following  sections. 

Figure  3(a)  also  shows  fhe  effecl  of  unload-reload 
cycling  and  hysteresis*  demonslraling  lhal  Ihe  material 
recovers  only  slighlly  after  Ihe  lirsl  cycle,  while  sub- 
sequenl  cycling  resulls  in  a  subslanlially  higher  relative 
recovery.  This  renders  Ihe  firsl  cycle  dislinclly  differenl 
from  all  subsequenl  loading  cycles,  as  has  been  widely 
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observed  for  CNT  foams. The  inset  in  Fig.  3(a)  shows 
the  zoomed-in  load-unload  segments  during  cycling  and 
indicates  that  both  the  amount  of  hysteresis,  corresponding 
to  the  energy  dissipated  per  cycle,  and  the  recovered 
displacement  after  unloading  diminish  with  each  sub¬ 
sequent  cycle.  Although  each  cycle  was  specified  to  be 
unloaded  to  only  10%  of  the  maximum  load  from  the 
previous  cycle,  in  some  cases,  the  applied  loads  were 
measured  as  slightly  negative,  likely  due  to  the  adhesion 
between  the  VACNT  tips  and  the  diamond  indenter.^^’^^ 

Several  indentation  load  (and  stress)-displacement 
curves  for  three  different  displacement  rates  are  shown 
in  Fig.  3(b).  Although  the  general  behavior  across  the 
different  rates  is  similar  in  appearance,  several  important 
distinctions  can  be  made.  For  example,  at  the  faster 
loading  rates  of  100  and  1000  nm/s,  there  appears  to  be 
a  rapid  extensive  displacement  burst  of  —20  pm  or  larger 
immediately  upon  attaining  ains,  resulting  in  a  temporary 
loss  of  contact  between  the  sample  and  the  indenter  tip. 
The  burst  size  appears  to  increase  with  loading  rate.  This  is 
in  contrast  to  the  gradual  load  removal,  where  contact  is 
maintained  at  all  times,  shown  in  the  tests  conducted  at  the 
slowest  rate  of  10  nm/s.  Based  on  the  in  situ  indentation 
tests  described  in  Fig.  2,  it  appears  that  the  burst  may  be 
caused  by  a  rapid  vertical  shearing  of  the  material  directly 
underneath  the  indenter  tip  along  the  indenter  edges. 
Indeed,  when  the  test  was  stopped  at  the  moment  of  the 
burst  [top  SEM  image  in  Fig.  3(b)  at  100  nm/s],  we 
observe  the  formation  of  a  large  imprint  of  depth  —19  pm 
mirroring  the  shape  of  the  indenter,  with  the  residual  depth 
equivalent  to  the  displacement  burst.  Closer  inspection  of 
the  imprints  indicates  that  the  material  directly  underneath 
the  indenter  is  mostly  unperturbed,  with  all  deformation — 
again — appearing  to  have  been  localized  in  the  vertical 
shear  region  along  the  rim  of  contact  between  the  sample 
and  the  indenter  tip. 

Further  loading  beyond  the  first  burst  is  characterized 
by  the  reappearance  of  the  sloped  plateau.  Such  sloped 
plateau  regions  have  been  observed  previously^’"  and 
were  attributed  to  an  inherent  property  gradient  in  the 
VACNT  microstructure  along  its  height.  Current  authors 
have  reported  a  similar  gradient  to  be  present  in  APCVD 
VACNT  films  where  bofh  fhe  CNT  number  and  mass 
densities  were  found  to  vary  as  a  function  of  the  film 
heighf.  If  is,  therefore,  reasonable  to  assume  that  a 
similar  density  variation  is  present  in  the  samples  de¬ 
scribed  here,  resulting  in  the  indenter  punch  encountering 
a  denser,  and  hence  stiffer,  material  in  the  deeper  regions. 
This  requires  the  application  of  a  higher  stress  to  form  each 
successive  buckle,  leading  to  an  increasing  global  slope 
in  the  plateau  region.  While  the  slopes  of  the  plateau 
region  appear  to  be  similar  among  all  three  rates  shown  in 
Fig.  3(b),  samples  deformed  at  the  faster  rates  sustain 
higher  stresses,  a  behavior  considered  typical  for  visco¬ 
elastic  solids. 


The  plateau  region  is  composed  of  several  undulations 
or  kinks  for  all  three  loading  rates.  Such  features  have 
been  observed  previously  during  in  situ  compression 
experiments  of  60- pm- tall  VACNT  bundles"  and  in  ex 
situ  tests,^’"’'^’^^’^"^’^^  as  well  as  more  recently  by  the 
present  authors  for  in  situ  compressions  of  30  x  30  pm 
height  X  diameter  VACNT  micropillars  made  from  the 
same  APCVD  process.^  The  in  situ  compression  studies 
of  Hutchens  et  al. ' '  revealed  that  the  formation  of  each 
undulation  corresponded  to  the  formation  of  localized 
folds  or  buckles  along  the  pillar  height,  with  the  faster  rates 
showing  drops  of  larger  amplitude.^’"  While — unlike  in 
compression  experiments — direct  observation  of  buckle 
formation  is  not  possible  during  indentation  in  the  interior 
of  the  sample,  it  is  reasonable  to  assume  that  a  simi¬ 
lar  deformation  mechanism  operates  during  flat-punch 
indentations.  Thus,  the  series  of  load  drops  in  the  load- 
displacement  curves  likely  correspond  to  the  coordinated 
sequential  lateral  collapse  events  somewhere  below  the 
indenter,  with  the  force  required  to  cause  them  being  much 
smaller  than  that  of  the  initial  instability. 

Unloading  from  a  depth  of  —65-70  pm  results  in  very 
little  recovery  as  shown  by  the  residual  imprint  in  the 
bottom  SEM  image  of  Eig.  3(b)  inset.  The  imprint  also 
shows  the  remnant  of  the  first  large  burst,  which  is  seen  as 
a  shear  line  at  a  depth  of  —20  pm  along  the  walls  of  the 
crater.  The  vertical  walls  of  the  indent  imprint  indicate  a 
clear  shearing  of  the  VACNTs  along  the  rim  of  the 
indenter. 

D.  Effect  of  substrate  on  buckle  formation 

Eigure  3  and  video  S 1  indicate  that  the  buckle  formation 
in  the  VACNT  film  is  strongly  influenced  by  the  con¬ 
straints  imposed  by  the  stiff  substrate.  In  the  taller 

—  142-pm-thick  sample  shown  in  Eigs.  1  and  2,  the  location 
of  the  initial  buckle  is  always  close  to  the  substrate,  while 
the  region  underneath  the  indenter  tip  remains  largely 
unaffected.  Eurther  loading  is  thought  to  cause  additional 
localized  folds  to  continuously  form  along  the  VACNT 
film  height,  and  the  indenter  starts  probing  the  folded/ 
collapsed  densified  region — an  effect  similar  to  that 
observed  during  compression  experiments  on  these 
VACNTs  patterned  into  micrometer-sized  bundles.®  It 
stands  to  reason  that  such  a  response  can  be  observed 
under  indentation  by  using  a  shorter  sample.  This  is  shown 
in  Eig.  4(a)  where  the  indentation  response  of  a  shorter 
25 -pm- thick  VACNT  film  is  compared  to  a  thicker 

—  142-pm-thick  film.  Both  tests  shown  in  Eig.  4(a)  were 
conducted  in  the  samples’  interior.  It  is  evident  that  while 
the  initial  loading  portions  of  both  samples  are  similar,  the 
shorter  VACNT  film  has  a  significantly  higher  slope  in  its 
plateau  up  to  an  indentation  depth  of  —19  pm,  suggesting 
a  higher  degree  of  stiffening  in  the  shorter  sample. 
Moreover,  the  plateau  region  of  the  shorter  sample  shows 
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FIG.  4.  (a)  Comparative  load-displacement  curves  at  100  nm/s  loading  rate  of  indentation  tests  conducted  on  a  thin  (25  pm)  VACNT  film  versus 
a  thick  (141.5  pm)  VACNT  film.  These  indentations  were  performed  on  the  interior  of  the  sample  far  away  from  edge  effects  using  a  120-pm-diameter 
cylindrical  flat-punch  indenter.  The  film  morphology  (b)  before  and  (c)  after  in  situ  SEM  indentations  using  a  60  x  80  pm  rectangular  flat-punch 
indenter  on  the  edge  of  the  sample  help  reveal  the  (c  inset)  progressive  buckling  seen  in  the  shorter  (25  pm)  VACNT  film  sample.  All  SEM  images 
except  the  bottom  (c  inset)  were  taken  at  86°  tilt  angle.  The  (c  inset)  image  was  at  60°  tilt  angle. 


more  prominent  undulations  than  during  indentation  of  the 
thicker  sample  at  the  same  rate.  Beyond  the  indentation 
depth  of  19  pm,  the  effect  of  the  stiff  Si  substrate  becomes 
significant  in  the  shorter  sample  as  displayed  by  a  steep 
rise  in  the  load-displacement  response.  Consistent  with 
these  observations,  stiffnesses  calculated  from  the  unload¬ 
ing  portions  of  the  load-displacement  curves  were  an 
order  of  magnitude  higher  than  those  of  the  thicker  film, 
showing  fhe  combined  effects  of  densification  and  of  the 
substrate  influence  in  the  shorter  Aims. 

To  gain  insight  into  the  possible  distinctions  in  mor¬ 
phological  features  between  samples  with  different  thick¬ 
nesses,  we  conducted  in  situ  indentation  experiments  in 
SEMentor  of  the  shorter  25-pm-thick  VACNT  Aim  sam¬ 
ple.  As  shown  in  Figs.  4(b)  and  4(c)  and  supporting  online 
video  S2,  these  tests  were  conducted  on  the  edge  of  the 
VACNT  Aim  [unlike  the  in-bulk  tests  shown  in  Fig.  4(a)]. 
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As  in  the  case  of  the  thicker  samples  (Fig.  2),  the  deformation 
extent  in  the  shorter  sample  is  localized  to  a  zone  directly 
underneath  the  indenter  tip  [compare  Figs.  2  and  4(c)]. 
However,  unlike  the  thicker  samples,  clear  evidence  of 
densiflcation  can  be  seen  in  the  shorter  VACNT  Aim  sample 
[Fig.  4(c)  inset].  Figure  4(c)  inset  also  shows  that  substrate 
proximity  has  resulted  in  multiple  folds  or  buckles  in  the 
shorter  VACNT  Aim. 

E.  Analysis  of  indentation  response 

Figure  5  shows  the  results  of  indentation  experiments 
conducted  on  the  ~142-pm-thick  VACNT  Aim  sample. 
Only  the  data  from  the  sample  interior  (i.e.,  no  on-edge 
tests)  are  shown.  Unloading  modulus,  deflned  in  Eq.  (1),  is 
plotted  as  a  function  of  loading  cycle  number  in  Fig.  5(a). 
This  plot  reveals  that  the  modulus  varies  within  a  range  of 
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FIG.  5.  (a)  The  unloading  modulus  (average  ±  standard  deviation) 
at  three  different  loading  rates  as  a  function  of  the  load-unload  cycle 
number.  The  lighter  shade  of  the  modulus  values  calculated  from  the 
first  unload  cycle  reflects  the  uncertainty  in  their  accuracy  due  to  the 
approximations  inherent  in  their  calculations,  (b)  %  recovery  as  a 
function  of  maximum  displacement  before  unload.  Unloading  the 
indenter  before  the  displacement  burst  (open  symbols)  shows 
significantly  higher  values  of  %  recovery  than  after  the  burst  (filled 
symbols),  (c)  The  loss  coefficient  values  as  a  function  of  load-unload 
cycles  showing  very  poor  energy  dissipation  under  indentation. 


120-175  MPa  as  a  function  of  displacement  rate  and  load- 
unload  cycle  number.  These  moduli  further  demonstrate 
a  strong  rate  effect,  with  the  fastest  rate  of  1000  nm/s 
consistently  resulting  in  the  highest  moduli  for  all  load- 
unload  cycles.  Such  a  rate  dependence  of  the  material 
stiffness  is  a  known  characteristic  of  viscoelastic  solids.^^ 
In  contrast,  repeated  cycling  appears  to  have  only  a  mar¬ 
ginal  effect  on  the  sample  modulus,  remaining  nearly 
constant  during  subsequent  cycles  for  each  displacement 
rate,  although  there  is  a  small  increase  from  the  first  to  the 
second  cycle.  We  note  that  the  above  method  of  calcula¬ 
tion  of  unloading  moduli  inherently  assumes  indentation 
of  an  elastic  half-space,  usually  represented  well  by 
a  monolithic  material.  Therefore,  they  may  not  be 
suitable  for  VACNT  indentation — ^particularly  for  the  first 
load-unload  cycle  where  a  block  of  material  separates 
from  the  outer  matrix.  In  spite  of  these  aspects.  Fig.  5(a) 
helps  to  underscore  the  higher  degree  densification  of  the 
VACNT  film  at  faster  rates.  We  also  find  that  the  higher 
modulus  at  faster  rates  leads  to  greater  instability  stress 
(Fig.  3),  which  increases  from  1.49  ±  0.2  MPa  at  10  nm/s 
to  1.54  ±  0.2  MPa  at  100  nm/s  to  1.75  ±  0.3  MPa  at 
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1000  nm/s  (average  ±  standard  deviation  over  10  or  more 
measurements).  Importantly,  the  indentation  moduli 
shown  in  Fig.  5(a)  match  well  with  those  measured 
during  compression  of  APCVD  VACNT  micropillar 
samples  [160-180  MPa,^  see  table  in  Fig.  7(b)]. 

Figure  5(b)  shows  the  trend  in  the  %  recovery  as 
a  function  of  maximum  displacement  before  final  unload 
across  the  three  loading  rates.  Two  distinct  regions  can  be 
identified  in  this  figure:  (i)  a  segment  of  higher  recovery 
between  50-80%  when  the  indentation  tests  are  stopped 
before  the  large  initial  displacement  burst  [open  symbols 
in  Fig.  5(b)],  while  (ii)  loading  beyond  the  burst  (filled 
symbols)  results  in  a  very  low  (<10%)  recovery.  This 
indicates  that  the  shear  event,  manifested  by  the  large 
displacement  burst  in  the  indentation  curve,  results  in 
permanent  damage  of  the  VACNT  structure  and  signifi¬ 
cantly  compromises  its  recoverability.  We  note  that  the 
virtually  nonexistent  material  recovery  from  postdisplace¬ 
ment  burst  during  indentation  is  in  stark  opposition  to  the 
almost  full  (—100%)  recovery  seen  in  these  materials 
under  compression  experiments^  [see  table  in  Fig.  7(b)]. 
We  did  not  observe  any  significant  deformation  rate 
effects  on  the  %  recovery. 

The  loss  coefficient  q  (see  Sec.  II)  is  shown  in  Fig.  5(c) 
as  a  function  of  load-unload  cycle  number.  The  data 
contain  only  those  results  for  the  tests  loaded  beyond  the 
displacement  burst.  As  a  consequence  of  the  marginal 
recovery  [Fig.  5(b)],  the  loss  coefficient  associated  with 
nanoindentation  into  this  material  is  also  very  low,  on  the 
order  of  0.003-0.006,  a  value  nearly  an  order  of  magnitude 
lower  than  those  for  compression  of  the  same  material^ 
[table  in  Fig.  7(b)]. 

The  viscoelastic  indentation  response  of  the  VACNT 
film,  in  terms  of  the  measured  values  of  their  storage 
modulus  (£”),  loss  modulus  (E"),  and  tanS  values,  is 
shown  in  Fig.  6.  Similar  to  the  %  recovery  results,  two 
distinct  responses  are  seen  in  this  figure  depending  on 
whether  the  indenter  is  oscillated  before  the  large  initial 
displacement  burst  (open  symbols  in  Fig.  6)  or  after  (filled 
symbols).  Both  E'  and  E'  increase  by  a  factor  of  —3^ 
once  the  displacement  burst  has  occurred  [Figs.  6(a) 
and  6(b)],  while  tanS  decreases  by  20-30%  after  the  burst 
[Fig.  6(c)].  Before  the  burst,  these  parameters  depend  more 
strongly  on  the  applied  load  during  oscillations,  as 
evidenced  by  the  larger  spread  in  the  data. 

We  find  the  storage  modulus  to  be  frequency  in¬ 
dependent  over  the  range  of  frequencies  probed.  On  the 
other  hand,  the  loss  modulus  (and  the  tanS)  strongly 
increases  with  increasing  frequency,  although  two  local 
minimas  exist  at  30  and  50  Hz  [Figs.  6(b)  and  6(c)]. 
Unfortunately,  the  cutoff  frequency  (50  Hz)  of  our  in¬ 
strument  prevents  further  study  of  this  behavior  at  higher 
frequencies. 

As  shown  earlier  (Figs.  2  and  3),  the  large  displacement 
burst  occurring  in  the  initial  loading  stages  results  in 
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significant  buckling.  Thus,  the  higher  E'  characteristic  of 
the  postburst  structure  indicates  that  the  material  is  capable 
of  storing  a  higher  amount  of  energy  in  this  state,  with  the 
excess  energy  presumably  being  stored  in  the  buckled 
regions.  The  corresponding  increase  in  E" — related  to  the 
dissipated  energy — is  at  a  lower  proportion,  thus  leading 
to  a  less  pronounced  viscous  response  (lower  values  of 
tan5)  after  the  burst.  The  post-burst  decrease  in  tan5  is 
likely  caused  by  the  increased  adhesion  between  the 
individual  nanotubes  in  the  buckled  regions  when  they 
are  brought  into  proximity  of  one  another  after  buckling. 
A  similar  response  has  been  observed  in  polymeric 
materials  like  rubber,"*^  where  an  increase  in  cross-linking 
leads  to  a  decrease  of  the  viscoelastic  response.  We  note 
that  the  trends  reported  in  Fig.  6  are  similar  to  those 
reported  previously  for  other  VACNT  microstructures 
under  indentation.*  The  storage  modulus  E'  shown  in 
Fig.  6(a)  matches  well  with  the  unloading  modulus  shown 
in  Fig.  5(a). 


0.14 
0.12 
I  0-1 

c  0.08 

(C 

0.06 
0.04 
0.02 
0 

,  .  0  10  20  30  40  50  60 

(C)  Frequency,  Hz 

FIG.  6.  The  viscoelastic  indentation  response  of  the  VACNT  film  at 
a  loading  rate  of  100  nm/s,  in  terms  of  the  measured  values  (average  ± 
standard  deviation)  of  their  (a)  storage  modulus,  (b)  loss  modulus,  and 
(c)  tanS  values,  as  a  function  of  frequency. 
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F.  Comparing  indentation  and  compression 

Figure  7(a)  displays  a  comparative  analysis  between 
indentations  into  the  ~  142-pm-thick  APCVD  VACNT 
films  in  this  study  and  the  compressions  on  30  x  30  pm 
(diameter  x  height)  APCVD  VACNT  micropillars  from 
our  previous  work.^  All  tests  were  conducted  at  the 
loading  rate  of  1000  nm/s.  The  indentation  stress  (aind) 
and  strain  (Sind)  values  were  calculated  using  Eq.  (3).  The 
following  observations  can  be  made  from  this  comparison: 
(i)  instability  stress  in  compression  is  —1.5  times  higher 
than  that  in  indentation,  (ii)  indentation  is  characterized  by 
a  significantly  larger  strain  burst  of  (Ac  =  0.3)  after  the 
instability  as  compared  to  compression  (Ac  =  0.05).  This 
can  be  explained  by  the  different  deformation  mechanisms 
associated  with  the  bursts.  During  indentation,  a  single 
global  shear  event  is  responsible  for  such  a  burst,  while 
under  compression,  a  relatively  small  burst  is  sometimes 
observed  just  after  yielding,  which  has  been  attributed 
to  a  temporary  loss  in  the  pillar-indenter  contact.^’" 
(iii)  After  yield,  the  slope  of  the  plateau  region  in  com¬ 
pression  is  significantly  higher  (by  a  factor  of  almost  6) 
than  during  indentation.  Thus,  while  the  compressive 
stresses  rise  from  —2.8  MPa  at  yield  to  4.8  MPa  at  the 
strain  of  0.6  (i.e.,  before  the  onset  of  densification),  the 
stresses  at  equivalent  strains  in  indentation  are  relatively 
constant,  in  the  2.0-2. 1  MPa  range.  Finally,  (iv)  after 
unloading  (to  10%  of  maximum  load),  the  material  under 
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FIG.  7.  (a)  Comparative  stress-strain  curves  of  VACNT  in  indentation 
and  compression,  (b)  Table  comparing  the  measured  values  from 
indentation  (this  study)  versus  compression  (Ref.  9). 
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indentation  shows  negligible  (<5%)  recovery,  while  com¬ 
pression  leads  to  an  almost  complete  (>95%)  recovery. 

While  determining  the  underlying  physics  of  the 
described  observations  requires  rigorous  multiscale  mod¬ 
eling,  some  of  these  phenomena  can  be  understood  in 
terms  of  the  differences  in  boundary  conditions  between 
the  two  loading  configurations  (indentation  versus  com¬ 
pression),  in  combination  with  the  overall  similarity  of 
the  VACNT  deformation  to  that  of  open-cell  foam 
materials."^®  The  characteristics  of  flat-punch  indentations 
into  VACNT  Aims  closely  resemble  those  reported  for 
metallic"^^’'*^  and  polymeric  foams. In  the  spirit  of 
foam-like  deformation,  the  indentation  zone  for  VACNTs 
is  confined  fo  fhe  region  direcfly  beneafh  fhe  indenfer, 
while  fhe  surrounding  regions  are  unaffected  (see  Fig.  2). 
The  lack  of  lateral  spread  in  fhe  deformed  region  is  fypi- 
cal  for  foams,  whereby  if  can  be  aflribufed  fo  fhe  non- 
exisfenf  Poisson’s  rafio.^’"*'^’"^^  This  is  diamefrically 
opposite  fo  monolifhic  materials,  where  the  plastic  zone 
underneath  the  indenter  tip  is  typically  represented  by 
a  continuous  hemispherical  plastic  zone.^^  Another  com¬ 
mon  characteristic  between  indentation  into  conventional 
foams  and  VACNTs  is  the  appearance  of  a  so-called  “tear 
line”^°’"^*  along  the  comer  of  the  flat-punch  indenter  where 
the  shear  force  is  largest.  The  depth  of  this  shear  line  in  a 
foam  generally  depends  on  its  strength."*'’’"^'  Unlike  foams, 
however,  the  constituents  of  the  VACNT  system — the 
individual  CNTs  themselves — are  nominally  vertically 
aligned  with  respect  to  the  indenter  tip  [Fig.  1(b)],  and 
thus,  the  vertical  plane  is  expected  to  be  the  plane  of  lowest 
shear  strength  in  the  VACNT  matrix.  Once  the  critical 
shear  stress  is  attained  during  loading,  the  shear-off 
proceeds  catastrophically  along  this  vertical  plane  through 
the  thickness  of  the  VACNT  film  down  fo  the  underlying 
substrate.  This  results  in  highly  localized  deformation 
along  the  shear  plane,  with  no  notable  perturbation  in  other 
regions.  Subsequent  loading  initiates  folding  or  buckling 
of  the  material  close  to  the  substrate  but  not  in  the  sheared- 
off  block.  In  the  extreme  case  of  indenting  a  short  VACNT 
film  (Fig.  4),  subsequenf  densiflcafion  of  fhe  already 
folded/buckled  CNTs  commences. 

On  the  contrary,  under  uniaxial  compression,  the  shear 
stresses  acting  on  the  VACNT  matrix  are  negligible. 
Instead,  the  applied  compressive  strain  is  accommodated 
entirely  via  the  formation  of  lateral  folds  or  buckles  along 
the  length  of  the  bundle,  while  the  remaining  portion 
remains  virtually  unscathed.^’"  Thus,  the  first  instability 
shown  in  the  compression  stress-strain  response  in  Fig.  7(a) 
is  related  to  the  onset  of  the  first  buckle  formation  in  the 
VACNT  matrix.  Further  loading  beyond  yield  results  in 
the  compression  of  the  collapsed  buckles.  The  resulting 
densiflcation  of  the  matrix,  in  combination  with  a  possible 
property  gradient  in  the  VACNT  structure,^’ “  results  in 
the  strong  positive  slope  of  the  compression  stress-strain 
curve  in  the  plateau  region. 


A  comparison  of  the  instability  stresses  between 
compression  (buckling,  2.69  ±0.1  MPa)  and  indentation 
(shear,  1.75  ±  0.3  MPa)  on  the  APCVD  VACNT  samples 
indicates  that  the  instability  stress  in  indentation  is  40% 
lower  than  that  of  compression  [Fig.  7(b)].  Interestingly, 
this  observation  is  opposite  to  that  of  traditional  indenta¬ 
tion  tests  into  monolithic  materials,  where  yield  in  in¬ 
dentation  typically  requires  1-3  times  higher  stresses  than 
for  the  uniaxial  case  for  most  solid^^  and  porous^^  as  well 
as  foam'^'’'^^’^®  materials.  The  two  distinct  modes  of  de¬ 
formation  in  VACNTs — shear  under  indentation  and  buck¬ 
ling  in  compression — are  likely  responsible  for  this  apparent 
reversal  in  their  instability  stress  values,  i.e.,  during  in¬ 
dentation,  the  APCVD  VACNT  film  reaches  ifs  critical 
shear  stress  before  it  can  buckle. 

While  our  findings  may  be  represenfafive  of  fhe  overall 
CNT  film  deformation,  if  is  imporfanf  fo  recognize  fhaf 
bofh  fhe  shear  sfrengfh  and  fhe  deformation  behavior  of 
the  VACNT  matrix  are  strongly  affected  by  its  micro- 
structural  details  such  as  tube  and  node  density  and 
tortuosity.  These  parameters  govern  the  entanglement  of 
the  individual  CNTs  within  the  matrix,  and  hence,  our 
observations  are  valid  only  for  the  particular  APCVD 
VACNT  film  sfudied  in  fhis  work.  Further,  even  the 
standard  APCVD  growth  techniques  for  VACNTs  are 
notoriously  difficult  to  control.  For  example,  despite  the 
growth  recipe  to  create  the  30- pm- tall  VACNT  pillars 
studied  in  Ref.  9  and  the  142-pm-tall  VACNT  Aims 
investigated  in  this  work  being  nominally  identical  (other 
than  the  growth  time),  it  is  reasonable  to  expect  substantial 
differences  in  the  VACNT  density,  alignment,  and  tortu¬ 
osity,  and,  especially,  their  respective  variations  along  the 
sample  heights.  Both  of  these  parameters  have  proved  to 
be  very  difficult  to  quantify  experimentally.  Similarly,  the 
two  VACNT  Aim  thicknesses  shown  in  Fig.  4  may  also 
show  a  similar  difference  since  they  were  grown  in 
different  batches.  Our  related  work  on  VACNTs  produced 
using  other  CVD  methods’*^  as  well  as  that  of  others  in 
the  field  have  yielded  a  wide  range  in  insfabilify  sfresses 
under  bofh  indenfalion  and  compression,  which  are  fhoughf 
fo  be  a  consequence  of  fhe  particular  details  (e.g.,  density, 
degree  of  entanglement,  etc.)  of  the  CNT  stmcture  and 
arrangement  in  their  matrix. 

We  attribute  the  large  difference  in  recoverability  of 
the  VACNT  matrix  in  indentation  (4.3%  ±  0.3%)  versus 
compression  (95.7%  ±  2.8  %)  to  the  different  deformation 
modes  and  boundary  conditions  under  these  two  loading 
modes.  The  almost  complete  recovery  of  the  VACNT 
structure  under  compression  in  Ref.  9  is  related  to  the 
elastic  unfolding  of  the  buckles  created  during  loading. 
This  has  been  explained  by  the  intertube  interactions  of  the 
CNTs  where  the  zipping  and  unzipping  of  CNTs  in  contact 
have  been  proposed  as  the  primary  criteria  for  the  superior 
viscoelasticity  and  creep  recovery  of  VACNTs. 

On  the  contrary,  recovery  in  indentation  is  inhibited  by  the 
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additional  constraints  of  the  surrounding  VACNT  material. 
Thus,  even  at  lower  strains  (before  the  large  shear  hurst), 
the  recovery  under  indentation  is  still  lower  [50-80%, 
Fig.  5(h)]  than  that  in  compression.  After  the  displace¬ 
ment  hurst,  recovery  of  the  VACNTs  in  indentation  is 
further  compromised  hy  the  permanent  nature  of  the 
damage  caused  hy  the  shearing  of  a  large  block  of  the 
material  from  the  matrix. 

Similarly,  the  viscoelastic  properties  of  the  VACNT 
matrix  are  also  different  in  indentation  (Fig.  6)  versus 
compression  (Ref.  9).  While  the  frequency  dependence 
of  tan5  is  similar  in  both  cases,  the  values  of  tanS  in 
indentation  (as  measured  after  burst)  were  20-30%  lower 
than  those  in  compression,  suggesting  a  lower  viscous 
response. 

IV.  CONCLUSIONS 

In  summary,  we  demonstrate  that  under  large  displace¬ 
ment  indentations,  APCVD-grown  VACNT  films  deform 
by  a  single  vertical  catastrophic  shear-off  event  along 
the  rim  of  the  indenter  tip,  manifested  by  an  exten¬ 
sive  ~20^0-pm-long  displacement  burst  in  the  load- 
displacement  curve.  Similarly,  to  foam  materials  with 
negligible  Poisson’s  ratios,  the  indentation  in  VACNTs  is 
also  exclusively  confined  to  regions  directly  beneath  the 
indenter.  Through  in  situ  nanomechanical  indentation 
experiments,  we  show  that  the  axial  strain  during  the 
catastrophic  shearing  event  is  accommodated  by  rolling  of 
several  coiled-up  CNT  sections,  or  “microrollers.”  We 
hypothesize  that  such  coiling  provides  a  low-friction  path 
for  deformation,  which  may  explain  the  lower  stresses 
associated  with  indentation.  Loading  beyond  the  initial 
shear  event  causes  a  20-30%  decrease  in  the  material’s 
viscous  response.  We  also  find  progressive  buckling  in  the 
VACNT  structure,  with  the  buckles  forming  primarily 
near  the  stiff  underlying  substrate.  Comparison  of  in¬ 
dentation  versus  compression  experiments  indicates  that 
APCVD-grown  VACNT  films  reach  their  critical  shear 
stress  before  buckling  in  indentation,  which  results  in  the 
catastrophic  shear  formation  in  their  microstructure.  As 
a  result,  these  materials  demonstrate  negligible  recovery 
in  indentation,  as  compared  to  their  almost  perfect  re¬ 
coverability  under  compression. 
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